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Using a new impurity density matrix renormalization group scheme, we establish a reliable 
picture of how the low lying energy levels of a S = 1 Heisenberg antiferromagnetic chain change 
quantitatively upon bond doping. A new impurity state gradually occurs in the Haldane gap as 
J' < J, while it appears only if J'/J > 7 C with I/70 = 0.708 as J' > J. The system is non- 
perturbative as 1 < J'/J < 7 C . This explains the appearance of a new state in the Haldane gap in 
a recent experiment on Y2- aJ Ca a ,BaNi05 [J.F. DiTusa, et al., Phys. Rev. Lett. 73 1857(1994)]. 
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Since Haldane [1] predicted that integer-spin antiferro- 
magnetic (AF) chains have a gap in the excitation spec- 
trum, much work [2,3] has been done on the 5 = 1 
^-^Heisenberg AF chain. In particular, the valence bond 
^^solid (VBS) model proposed by Affleck et al [3] has suc- 
OVeeded in verifying this prediction. Recently, very accu- 
^rate estimates of the Haldane gap (A = 0.41049(2)J) 
bJhave been achieved by White and Huse [4] and Go- 
^linelli, et al [5], using the density matrix renormaliza- 
•^T^tion group (DMRG) method and exact diagonalization 
O^together with proper extrapolations, respectively. Evi- 
^v^dence for the gap has also been clearly shown by some 
experiments on quasi- ID 5 = 1 AF materials: NENP [6] 
^vqand Y 2 BaNi0 5 [7-9]. Impurity effects have not been fully 
(^studied yet, although the VBS picture indicates that the 
i-Hphysics of low lying excitations could be dramatically af- 
OQfected by the introduction of impurities. ESR on NENP 
^^[10] shows evidence for the existence of free 5=1/2 
Q^spins at the chain ends upon magnetic doping with Cu 
"^or non-magnetic doping with Zn, Cd and Hg, whereas the 
^specific heat measurements on Y2BaNi 1 _ a; Zn a ;05 show a 
£nSchottky anomaly consistent with free 5=1 spins [9]. 
I Recently DiTusa et al [8] have performed an experiment 
^with carrier doped Y 2 _ :E Ca :E BaNi05. The inelastic neu- 
Qtron scattering shows that the Ca doping produces new 
^states in the Haldane gap. According to polarized x-ray 
absorption spectra, the added hole is localized on the 
oxygen site with 5 = 1/2 spin, which might either di- 
rectly break the exchange coupling J and thus couple to 
5 = 1 Ni +2 ions or otherwise change J into J 1 . It is clear 
that the impurity effects are fundamental for further un- 
derstanding the properties of Haldane gap systems. 

Recently, S0rensen and Affleck [10] have studied the 
low-lying excitations in these doped Haldane systems us- 
ing the DMRG method [11] and a free massive field the- 
ory. One of the most surprising findings is that critical 
values for J' are required to induce impurity states in 
the Haldane gap both in the weak and strong coupling 
regimes. However, as remarked by the authors, their cal- 
culation using the standard DMRG method becomes pro- 



gressively inaccurate for small or large couplings, since 
the system is split exactly at the impurity bond and cor- 
relation effects of the impurity bond entering truncated 
blocks are enhanced as the chain grows. Their results do 
also not agree with those obtained by a Schwinger boson 
technique [13], which gives rise to an impurity state in 
the middle of the Haldane gap being almost independent 
of J'. Therefore, a substantially more accurate calcula- 
tion is needed to clarify the fundamental properties of 
such doped Haldane systems. 

The DMRG technique [12] invented by White has suc- 
cessfully provided highly accurate results for some homo- 
geneous systems. For a given number of states kept in 
the truncation procedure, numerical results have a higher 
precision using open rather than periodic boundary con- 
ditions. On the other hand, restricting oneself to open 
boundary conditions leads to an unfavourable scheme for 
the application of DMRG to impurity systems. One can 
nevertheless provide very accurate results even for impu- 
rity systems if one is able to work with a periodic bound- 
ary and sufficiently many states. Periodic boundary con- 
ditions allow for a proper placement of the impurities 
in the chain. A general scheme suitable for an impurity 
problem (called the impurity scheme hereafter) is shown 




FIG. 1. Impurity DMRG scheme 

in Fig. 1. In this scheme, a truncation is carried out only 
for the blocks B so that the impurity part contained in 
the dashed box is not subject to the truncation proce- 
dure [15]. Moreover, with increasing chain length, added 
sites depart further and further from the impurity part, 
which is in the spirit of Wilson's renormalization group 
as applied to the Kondo impurity problem in momentum 
space. In order to keep sufficiently many states, a well- 
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optimized code is designed, which is suitable for quite 
general DMRG calculations. 

In this letter, we consider the case where due to dop- 
ing one bond J of a closed 5 = 1 Heisenberg AF chain is 
changed to J'. The corresponding Hamiltonian is 

L-l 

H= J ^ S, • S, +1 + J'Sx • S L (1) 

where J > and < J' < oo. It is essential to 
point out some relevant properties of eq. (1) for special 
cases in order to understand our numerical results: i)For 
7 = J 1 1 J = 1, eq.(l) reduces to a 5 = 1 Heisenberg AF 
chain of length L [2]. For even L, the ground state is a 
singlet. In the thermodynamic limit continuous excita- 
tions start at the Haldane gap and the lowest one is a 
triplet state with 5 = 1 and q = 7r. ii)As 7 = 0, i.e. the 
open boundary case [14], the lowest excitation with 5 = 1 
becomes degenerate with the singlet ground state in the 
thermodynamic limit. The Haldane gap still exists, but 
is now the energy difference between the ground states 
and the 5 = 2 states. In this case, the chain end spins 
are fully polarized, iii) As 7 — » 00, the spins linked by the 
impurity bond form a singlet, which leaves the problem 
of a L — 2 sites chain with an open boundary. All proper- 
ties of the open chain are present apart from this isolated 
singlet. As long as 7 7^ 0, 1, 00, we believe that the lowest 
excitations are rearranged and an impurity state is asso- 
ciated with the lowest excited state with 5 = 1. In the 
following, we employ exact diagonalization and DMRG 
methods together with proper extrapolations to clarify 
how the triplet state descends from the bottom of the 
continuum spectrum at 7 = 1 to the ground state at 
7 = 0, 00 and whether the properties of the system in 
the case of 7 > 1 are the same as for 7 < 1. 

Now we turn to some important details of our nu- 
merical calculation for impurity energy and more details 
will be given elsewhere. We use the standard DMRG as 
7 = 0, 1, for which 500 and 1100 states are kept in each 
block B, respectively. In the impurity scheme, we con- 
duct the calculation exactly up to the 16 site chain which 
contains 729 states in each block B, and subsequently we 
keep 600-630 states for the remaining calculation. For 
7 = 1, we found that the impurity scheme keeping 600 
states gives a result of similar accuracy as that given 
by the standard DMRG scheme using periodic boundary 
and keeping 1100 states. As 7 7^ 0,1, we also checked 
how reliable the standard DMRG results are using either 
open or periodic boundary. We found that the inaccuracy 
becomes apparent as 7 is either small or large, whereas 
the impurity scheme gives very accurate results even for 
very small and very large values of 7. Up to 40 site- 
chains [16], the truncation error is at most 7.0 x 10~ 9 
and typically of the order of 10~ 10 . 

To determine the accuracy of our calculation, we 
checked it by keeping different numbers of states. The 



systematic error of the energies due to truncation is 
smaller or of the order of 10 -6 . Moreover, exact re- 
sults for 7 = 1 are available for up to 22 sites, accord- 
ing to which the difference to DMRG results with 1100 
states kept are 4.0 x 10~ 9 and 4.3 x 10~ 8 for ground 
state energy and Haldane gap, respectively. We obtained 
A = 0.410499 from DMRG of 22 sites with Shanks tran- 
formation, which is essentially the same as that given by 
exact diagonalization with Shanks transform. [5] To ob- 
tain final results we eliminate finite size effects for each 
7 by both fitting all data to exponential laws and using 
Shanks transform with smaller chains (L = 22 ~ 26). We 
found that these two results coincide very well. The er- 
rors of the final results we obtained for the impurity state 
are at most in the fourth digit as is partially shown in 
Table I, which we consider sufficiently accurate to ensure 
the reliability of our conclusions. 
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FIG. 2. Energy versus J' /J (left) and J/ J' (right) 

The dependence of the energy of the low-lying states 
on bond doping is given in Fig. 2. The horizontal axis 
indicates the ground-state energy. The darker shaded 
area represents the continuum spectrum and the straight 
line at its bottom marks the Haldane gap which is inde- 
pendent of the impurity coupling J'. According to the 
analysis of finite chains of up to 16 sites, we found that 
the bottom state has 5 = 2 at 7 = 0, 00 and 5 = 1 oth- 
erwise. Using DMRG, we checked A = 0.411 for 7 = 0.6 
[17]. Solid curves represent the impurity state, whose ap- 
pearance in the Haldane gap depends on the strength of 
the impurity coupling 7 J. In the weak coupling regime 
shown on the left side, it appears as long as < 7 < 1, 
whereas in the strong coupling regime shown on the right 
side it does only if 7 C < 7 < 00. 

To interpret the above results, we show some detailed 
data in Figs. 3 and 4 and insets show amplified critical 
regimes. For any finite L, the impurity state has 5 = 1 
forming a degenerate triplet. In order to characterize this 
state, we calculated the energy difference A i (7) between 
the ground state and the triplet state. The energy of this 
triplet state can be evaluated as the lowest energy in the 
S z = 1 subspace. In Fig. 3 A i (7) is denoted by an open 
square for given L and 7. The values of A L (7) decrease 
from top to bottom as L increases for given 7. For each 
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L G [8, 22] a fitting curve links these squares as a guide 
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FIG. 3. Weak coupling region 

to the eye and dashed lines denote A i (l). As is explicitly 
shown in Fig. 3, we find 

A L ( 7 < 1) < A L (1). 

This behaviour certainly persists as L — » oo. Using 
the extrapolation discussed above, we obtain A°°(7) as 
denoted by solid dots, which is fitted to a solid curve. 
This curve is below the Haldane gap denoted by a solid 
straight line. However, the situation is completely differ- 
ent for 7 > 1 as shown in Fig. 4. For any finite L, 




J / J" 
FIG. 4. Strong coupling region 

there exists an overshooting and the fitting curve inter- 
sects twice the dashed line which indicates the Haldane 
gap at that particular L. It is clearly seen that 

A L ( 7c L < 7) < A L ( 7 = 1) < A L (1 < 7 < 7 C L ), 

where 7^ is the value of 7 corresponding to the inter- 
section point for given L. To estimate 7 C for the critical 
point, we first utilize an extrapolation of the intersection 
points 7^. Secondly we determine the point at which 
the fitted curve of A 00 (7) drops below the Haldane gap. 



These two independent estimates yield 1/7 C = 0.708(2) 
where the error shows their difference. 

In order to further investigate the properties of such a 
doped system, one can map the original problem onto a 
L — 2 system by projecting out the singlet spins as 7 ^> 1 
[11]. The leading term in I/7 is still an impurity bond of 
a strength 4 J/37 and other terms are 0(1/ j 2 ). Such a 
mapping certainly breaks down as 7 < 4/3, where either 
higher orders or the break up of the singlet become cru- 
cial, resulting in significantly different properties for the 
strong coupling as compared to the weak coupling regime. 
Therefore, according to our numerical results, the prop- 
erties of such a doped system are non-perturbative as 
1 < 7 < 7c ) otherwise the problem can be treated prop- 
erly in 7 or I/7 in the weak or strong coupling regime, 
respectively. 
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FIG. 5. Spin polarization Pi 



Furthermore, it is instructive to analyze the spin po- 
larization induced by J' 7^ J. In this case, we calculate 
Pi = < ll^^ggll >, i.e. the expectation values of Sf +35 
in the state with S z = 1, and the impurity bond is trans- 
lated to the center of the periodic chain. As 7 = 1, 
we have P; = 1/L. In Fig. 5, we show Pi for four 
values of 7 with L = 70 and 300 states kept. In the 
weak and strong coupling limits, the Pi are essentially 
the same apart from a singlet in the middle as shown in 
(a) and (b) for 7 = 0.1,40/3, respectively. As 7 = 1/0.8 
which is in the non-perturbative regime, the region of 
sizeable polarization (c) is still similar to the cases (a) 
and (b). On the other hand, we find for (a) and (b) 
x(7,«) = lim l _> oo L P(i) very small for i ^> 1 while in 
the critical regime (c) this limit gives 1. In this sense 
we can distinguish localised from delocalised states. As 
can be seen from (d) for 7 = 1/0.6, x(7,i ^> 1) is again 
smaller than 1 but the change is smooth. We believe that 
x(7, i 1) behaves analogous to A°°(7), which may give 
an independent way to determine the critical point. 

Now we discuss the properties of such a doped system 
as 7 — » 7^ and 7 S with 7^ = 0, 00 and 7 S = 1~,7+, re- 
spectively. Based on the above analysis and perturbation 
theory for 7 — » 7u) , one finds 
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A oo W _ / a w J~ ( + 0(~i 2 ) 
17J " \ a s J/ 7 + 0(l/ 7 2 ) 



7< 1 
7 > 1. 



We extracted a„, = 1.1337 from 7 = 0.025, 0.05 and a s = 
1.5229 from 7 = 20, 40, which gives a s /a w = 1.3433. The 
deviation from 4/3 reflects higher order contributions in 7 
or 1 /'j. On the other hand, according to the VBS picture, 
a free S = 1/2 edge state exists as 7 = 0. One can project 
the two spins linked by the impurity bond onto a S = 1/2 
respresentation [11], giving A°°(7 <C 1) = a 2 J7 with 
< \\Sl |1 >=< 1|5£|1 >= a/2 as 7 < 1. Consequently 
one has a 2 = a w , which gives rise to < l|Sf |1 >= 0.5324 
being essentially the same as that obtained by previous 
calculations under open boundary conditions [12,18]. Re- 
garding the properties of the system near 7 S , one may as- 
sume a power law such as A°°(7) = A + a n J{^ s —7)™ as 
7 — » 7 S [11]. According to the perturbation theory with 
respect to J' — J for 7 S = 1 _ , we can exclude n = 1 since 
«i = A i (7 = 1)/L. By subtracting this term, we obtain 
a-2 = 0.0417 with n = 2 using the results for 7 = 0.95 
and L = 40. Fitting the data to a polynomial in 1/L, 
we have a? ~ 0.000(4) being zero within the error of the 
calculation, so that n > 3. Clearly both more analytical 
and numerical work are needed in order to completely 
determine this critical behavior. 

In conclusion, we have quantitatively established a re- 
liable picture of the appearence of the impurity state in 
the bond doped Haldane system concerning the off-chain 
doping material Y2- x Ca a; BaNi05, where a new state oc- 
curs in the Haldane gap. In distinction to Ref. [11,13], 
we found that in the weak coupling regime, an infinitely 
small deviation of J' from J gives rise to such a triplet 
state, whereas in the strong coupling regime, the system 
exhibits as a fundamental difference not only the forma- 
tion of a singlet for the spins linked by the impurity bond, 
but also the existence of a critical point. As 1 < 7 < 7 C , 
the system is clearly non-perturbative. 
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helpful discussions. 
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TABLE I. Singlet-Triplet energy difference as function of 7 
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